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ABSTRACT 


An "absorbing-sphere" model based on the Landau-Zener method was 
developed for calculating the energy dependence of the total cross sections 
for the reaction 

A + + b” A* + B + AE 

in which A or B may be atomic or molecular species. It was found that 
the cross sections were highly dependent on the electron detachment 
energy of B . In a second research task, the ionization (Penning and 
associative) total cross sections for Ne (2p 5 3s 3 P) 4* Ar were calculated 
over a large energy range and compared with the experiments of Micha, 

Tang, and Muschlitz. 
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INTRODUCTION 


The main objectives of the research under this contract are to 
perform fundamental studies in the field of low-energy electronic, ionic, 
atomic, and molecular impact phenomena. The theoretical studies are 
related to this laboratory’s experimental investigations of low-energy 
collisions and processes that are important to the understanding of the 
upper atmosphere. This research is a continuation of the work under 
Contract No. NASr-49(07) ; SRI Project No. -4070. 

SUMMARY OF RESULTS 

During the period of this contract, two separate research tasks 
were completed. The first, which was also partially supported by a DNA 
project, consisted of theoretically predicting the cross section energy 
dependence for ion-ion recombination of molecular species. The results 
were compared with existing experimental measurements. A second research 
task consists of trying to understand the unexpected energy dependence 
of the ionization cross sections measured by Micha, Tang, and Muschlitz 
on the Ne* + Ar system. 

Both of these works have just recently been submitted for publication. 
A summary of the results is given below, and a more detailed description 
is given in the appendices. 

A. Ion-Ion Recombination Total Cross Sections 

A procedure has been set up for calculating the total cross 
sections for the ion-ion recombination reaction 

4 * — * 

A + B -* A + B + AE • 
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The method was developed especially for the ca nes where A or B 
are complicated molecular species and where more exact calculations 
are not possible* The procedure is based on the Landau-Zener 
method for calculating the inelastic cross sections, A semiempirical 
method has been developed 1 with which to estimate the necessary 
coupling matrix elements as a function of the crossing distance 
R^ so that calculations can be made (the coupling matrix element 
research was supported by NASA by the contract that preceded this 
one) 0 We have previously used the Landau-Zener method to compute 

the ion-ion recombination cross sections as a function of energy 

+ — + —+ — + — 
for the systems H + H , He + H , 0 +0, and N + 0 , and we 

have found it to work surprisingly well 2 (this research was also 

supported by the preceding NASA contract). The multistate 

Landau-Zener calculations included many states ana, as an example, 

+ — 

in the N + 0 reaction, there were approximately 25 channels 
(excited electronic levels) available for reaction* If we 
proceed to molecular systems, we find that the number of available 
channels increases dramatically because of the large number of 
electronic, vibrational, and rotational states available. 

This leads t's to the assumption that for the more complicated 

systems we may approximate the physical situation by using a high 

density of crossing states. Then by using the Landau-Zener 

method and the semiempirical coupling matrix elements, we may 

obtain a critical crossing distance within vhicn there will 

be reaction. From R , the total cross sections can be calculated 

c 

by assuming an absorbing-sphere model, i.e., unit probability 

for reaction within R . 

c 
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For the systems studied experimentally, we find good agreement between 
theory and experiment for cases where there is an atomic negative ion, but 
for the molecular negative ion systems we find the theory tends to under- 
estimate the magnitude of the reaction rate obtained from extrapolation 
of the merged-beam results. In the reaction rate calculations, the cross 
sections are mainly determined by the electron detachment energy of the 
negative ion with lower electron detachment energies producing larger cross 
sections (the degree of excitation of the positive ion plays only a small 
part). We believe that since experimentally the negative ions are formed 
in a duo-plasmatron under very hot conditions, it is not unlikely that the 
molecular ions are formed in excited electronic or vibrational states. 

This would lead to an observed cross section that is larger than the 
calculated one, which is based on the assumption that the ions are in their 
ground states. This would also explain the reason why the more complicated 
atomic systems are in agreement with theory since the excited levels are 
not available to atomic negative ions. Further work, both experimental 
and theoretical, will be needed to verify these predictions. 

A complete description of the above work is presented in Appendix A 
and has been accepted for publication in The Journal of Chemical Physics 3 . 

B. Ne* + Ar Ionization Total Cross Sections 

The recent experimental data of Micha, Tang, and Muschlitz 4 on 
the low energy (0.01 to 0.1 eV) dependence of the ionization total cross 
sections for the sum of the Penning ionization process 

Ne ( 3 P 2 , 3 P 0 ) 4* Ar — * Ne + Ar + + e (l) 

and of the associative ionization process 

Ne ( 3 P 2 , 3 P 0 ) + Ar NeAr + + e ( 2 ) 
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display an unexpected energy dependence that has not been explained by 
theoretical calculations. These authors find cross sections that first 
decrease as the energy increases, then pass through a minimum, and later 
increase at the highest energies measured. Calculations by these authors 
using a Lennard- Jones, L-J, (l2,6) interaction potential and a short range 
coupling width r(R) were able to reproduce the low energy, monotonically 
decreasing behavior of the cross sections but did not follow the change 
in slope. Further experimental work by Tang, Marcus, and Muschlitz 5 on 
Ne + Kr and Ne + Xe show the same type of behavior. On another noble 
gas system — He (2 3 S and 2 i s) + Ar — the thermal energy cross section is 
comparable to or smaller than the high energy results at 100 to 1000 eV. 

It therefore appears that the shape of the energy dependence of the 

Ne + Ar cross sections could be common to many of the noble gas ionization 

reactions. 

* 

We have found that the Ne + Ar data can be reproduced with the 

same L-J (12,6) Potential that is used in Ref. 4 if the coupling width 

has a more realistic exponential form. Furthermore, calculations at high 

energies, E « 100 eV, indicate that the ionization cross section is of 
° 2 

the order of 10 A . This is similar in magnitude to the high energy 
experimental results for He + Ar that have been measured here 6 . 

Since the experimental measurements are most conveniently made at 
thermal energies, E <0.1 eV, and at high energies, E ^ 100 eV, the 
theoretical technique provides a means of interpolating the cross sections 
to the energy region 0.1 to 1.0 eV, the region of interest for upper 
atmosphere processes. 

A more complete description of this work is given in Appendix B. 

The manuscript has been submitted to Chemical Physics Letters for 
publication. 
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ABSTRACT 


An "absorbing-sphere" model based on the Landau-Zener method is 

o 

set up for calculating the upper limit thermal energy (300 K) reaction 
rate and the energy dependence of the total cross sections for the 
ion-ion recombination reaction 

+ - ^ 

A + B — * A + B + AE 

in which A or B may be atomic or molecular species. The crucial parameter 
needed for the calculation is the electron detachment energy for the outer 
electron on B . For a given system and energy, we find that, as the electron 
detachment energy is decreased because of the existence of metastable elec- 
tronic states or excited vibrational states, the cross sections increase. 

The calculated cross sections and rates are compared with afterglow results 
and with those obtained from "merged-beam" experimental data. We find good 
agreement with the existing results for atomic negative ions. There are 
differences for molecular negative ions, and the possible cause of these 
differences is discussed. 
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INTRODUCTION 


The theoretical development for the calculation of the total cross 
sections for the ion-ion recombination reaction 

A + B —►A +B + ( 1 ) 

may take one of two principal paths. The first, and most exact, is to 
calculate the ab initio potential energy surfaces for all the electronic 
states of reaction (l) and then to perform coupled-channel calculations 
to obtain the appropriate transition probabilities from which total 
cross sections may be computed. This method is time consuming, and the 
calculations have been performed only for the simplest reaction, that of H -f H 
The second approach is to use semiempirical methods to calculate the 
parameters needed to obtain reasonable estimates of the cross sections. 

This is the procedure that is described here. We will be concerned with 

o 

estimating the thermal energy (300 K) reaction rate and the energy depend- 
ence of the cross sections for reaction (l) when A or B or both may be 
atomic or molecular species. In particular, we will be concerned with 
species of upper atmosphere importance, since reaction (l) is an important 
mechanism in the final removal of charge in the D region of the atmosphere. 3 

The procedure to be described is based on the Landau-Zener method 
for calculating the inelastic cross sections. 4 “ 7 A semiempirical method 
has been developed 8 with which to estimate the necessary coupling matrix 
elements H 12 as a function of the crossing distance so that calculations 
can be made. We have previously used the Landau-Zener method to compute 
the ion- ion recombination cross sections as a function of energy for the 

+ — + •» -f* — _ 

systems H + H , He +H,0 +0, and N + 0 , and we have found it to 

work surprisingly well. 9 / 10 The multistate Landau-Zener calculations 
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included many states and, as an example, in the N +0 reaction, there were 
approximately 25 channels (excited electronic levels) available for reaction. 
If we proceed to molecular systems, we find that the number of available 
channels increases dramatically because of the large number of electronic, 
vibrational, and rotational states available. 

This leads us to the assumption that for the more complicated systems 
we may approximate the physical situation by using a high density of cross- 
ing states. Then by using the Landau-Zener method and the semiempirical 
coupling matrix elements, we may obtain a critical crossing distance 

within which there will be a reaction. From R , the total cross sections can 

c 

be calculated by assuming an absorbing-sphere model, i.e., unit probability 

for reaction within R . 

c 

We have compared our results with those obtained from an extrapolation 
of merged-beam data to thermal energies 9-14 and also with some afterglow 
measurements. 15 ; 16 In the systems studied experimentally, we find good 
agreement between theory and experiment for cases where there is an atomic 
negative ion, but for the molecular negative ion systems we find the theory 
tends to underestimate the magnitude of the reaction rate obtained from 
extrapolation of the merged-beam results. Possible reasons for this dis- 
crepancy are given in a later section. 

The energy dependence of the total cross sections may also be estimated 
by using the absorbing-sphere model. As in the reaction rate calculations, 
the cross sections are mainly determined by the electron detachment energy 
of the negative ion with lower electron detachment energies producing larger 
cross sections. Comparisons are made with experimental merged-beam data for 
several 0 systems, and the energy dependence of the cross sections is 
found to be well reproduced. The theoretical cross section energy dependence 
for several other negative ion systems is also given. 
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THEORY 


To determine the critical distance R within which there is 

c 

react ion^ it is best to review the application of the two-channel 
Landau-Zener method to reaction (l). It will then be possible to 
ascertain the functional dependence of the crossing parameters on the 
cross sections and to extend these results to the multicrossing case. 

In the two- state Landau-Zener approximation, the total cross 
section is given by (all quantities will be in atomic units) 

Q 12 = 4tiR x 2 (1 + AE/E) F S U) , (2) 

where F a (\) is the tabulated exponential integral 5 ; 6 

00 

F 3 (\) = J e~^ x (l - e * x )x 3 dx , (3) 

1 

and X is related to the velocity v and crossing parameters by 



In Eq. (4), H l2 is the coupling matrix element evaluated at the 

crossing distance R^, and V ; is the derivative of the potential with 

respect to R at R . For ion-ion recombination in which we have a 
x 

Coulomb potential crossing a dipole-induced dipole R“ 6 potential at 

X 4 rj 

large distances, R ~ 15 a.u., it is a good approximation to let 


A - 4 




and 

AE = -V . ( R ) = R _1 (6) 

xx x 

Substituting Eq. (5) into Eq. (4) ; we have 
, nJT’tt R 2 H 12 2 (R ) 

X = Jj x x (7) 

(E + AE) 1 / 2 

We may make a further approximation for the special case of thermal energy. 
Since the important range of crossings lies between 15 and 50 a 0 , we find 
from Eq. (6) that AE lies between 0.5 and 2.0 eV. Therefore, at thermal 
energies AE » E, so that by substituting Eq. (6) into Eq. (7) we arrive 
at the thermal energy result for X: 

X = s/2* TT R 5 / 2 ^/ 2 H 12 2 (R ) . (8) 

X X 

In the two-channel case, F 3 (X) has a maximum value of 0.113 when X 
is equal to 0.424. 

If we now extend these results to the multichannel case where 

there are a large number of states available, we may expect complete 

reaction for values of R less than R where R is the value of the 

c c 

crossing distance that satisfies X = 0.424 in Eq. (8). The total 
-cross section from Eq. (2) will then be 

Q = m c 2 [1 + (R c E) -1 ] . (9) 

To determine X, we must know its dependence on R. The coupling 
matrix element H 12 is the only unknown, but from previous work on 
one-electron transfer, we have been able to parametrize it in terms 
of the effective ionization potential of the reactants and products. 8 
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We may write the coupling matrix element in the reduced form 


H? 2 = 1.044 R* exp(-0. 857 R*) 


where 


H ? 2 = 




t 


( 10 ) 


and 


_ (a + y)R 

2 


( 11 ) 


The quantity a equals (21^) , where 1^ is the detachment energy of the 
electron being transferred, in a.u. For a ground state negative ion, 1^ is 
equal to the electron affinity. The quantity y equals , where I f is 

the effective ionization potential of the transferred electron in its 
product state. Equation (lo) is valid only at large distances, R > 10 a.u., 
where the electron is transferred through the potential barrier determined 
by the exponential tails of the wave functions for the reactant and product 
states and for values of ^ and 1^ greater than ~ 0. 4 eV. 


By knowing the electron detachment energy of the negative ion 
and by assuming that there is a high density of states available for 
reaction, we may then determine y as a function of R by using the 
relationship 


1 1 ml i + R “‘ (12) 

Using Eqs. (lo) and (ll), it is now an easy matter to obtain X, 

Eq. (7), as a function of R. From numerical calculations using 
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large numbers of final states, we find that the value of X to use 

for determining R for the multichannel case is approximately 0.15. The 
c 

value \ = 0.424, which is correct for the two-channel case, tends to under- 
estimate the value of the cross section by about 10$. The use of X = 0.15 

to determine R has been verified by Landau-Zener calculations employing 
c 

large numbers of product channels. For a realistic system where there 
obviously is not a continuum of product states available, the model 
continues to work since each state can now accommodate a higher percentage 
of the total reaction probability due to less competition with neighboring 
channels. The cross sections are reduced by only 20$ when there are as 
few as 5 product states. 

The formulas presented above are contingent on a high density of 
states being available for reaction. In reality, there may be significantly 
fewer, so that the cross sections calculated from the above formulas should 
be an upper limit on the true situation. 

Because of the Coulomb attraction, the product E*Q approaches a 
constant at thermal energies, Eq. (9), so that the reaction rate a , which 
is an average over a Maxwellian distribution of ion velocities, becomes 



For molecular systems, the coupling matrix elements, Eqs. (lo) and 
(ll), should be modified by the Franck -Condon factors q that represent 
the overlap of the vibrational energy levels to read 4 

H? 2 - 1.044 q x / 2 R*exp(-0.857 R*) . (l4) 
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The Franck-Condon factors are quantities less than unity, so again we find 
that we have an upper limit on the cross section. For the cases to be 
studied, even if all were set equal to 0.1, the cross sections would be 
decreased by only about 20$. 

Within the absorbing-sphere model described above, we find that the 
cross section is dependent only on the value of the electron detachment 

energy of the negative ion and the reduced mass of the system. A similar 
result is obtained in the work of Radtsig and Smirnov 18 who use a different 
approach to this problem. 

In a given system, as the electron detachment energy is decreased 
because of vibrational or electronic excitation, the reaction rate 
increases. Physically, this corresponds to the outer electron on the 
negative ion being more loosely bound, making it possible for transfer 
to occur at larger intermolecular separations. If the negative ion is 
present in an excited electronic or vibrational state, we would predict 
that the reaction rate will generally be greater than when the negative 
ion is in its ground state. This problem usually arises with molecular 
negative ions, but there may be some difficulty with an atomic negative 
ion system if the negative ion can be present in a metastable electronic 
state, such as C . 

To remove some of the objections to using the Landau-Zener model, 
namely the tunneling problem, we have compared the Landau-Zener transition 
probabilities with the close-coupled values tabulated in reduced form by 
Delos. 19 We find that the cross sections obtained by the Landau-Zener 
method are below the exact results and should be increased by approximately 
10$. This factor has been included in our results. 
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TI31HLIAL ENTERGY REACTION RATES 



he above formulas were first applied to systems in which the 
ve icn was ar: atonic state. For these cases, the electron detach- 
er gy is equal to the electron affinity. Since this quantity is 
r.cv:‘. : the calculations are straightforward. The results for several 
v are given in Table X, In Table I, the * limits on the theoretical 
on rates correspond to what we consider to be reasonable limits on 
curacy of the 11 2 generator. By comparing the matrix elements with 

rg experimental and theoretical data of H 12 vs R for one~elect?cn 

x 

sr cazs-i, v/o find that 35^ of the data lie within a factor of 3 
form given in Fqs 0 (ic) and (ll). The limits shown on Table I 
port, co 'she iacxor of 3 limits on H 12 . 


or the 


surpr 


atom-atom systems, we find good agreement with the extra- 

— + 

cjv: beam results. One exception, however, is the H ? K 
e, the theory overestimates the reaction rate. This result 
isiiig, however, since for this system at thermal energies, 
nnel 'the n~3 excited state of H) is open for reaction* The 


f ion that 
on is inva 


rate cc 


mere are a large number of states available for 

id in this case, leading us to obtain an upper limit 

siderahly above the experimental result. For the 
+ 


,0 4* X , and 0 + FT a reactions, the agreement with experiment 

cu Ir, these cases, we would expect reasonable agreement because 
are a large number of channels available for reaction. 
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TABLE I 


ION- ION THERMAL ENERGY RECOMBINATION REACTION RATES: 
ATOMIC NEGATIVE ION 


System 

Electron 
Detachment 
Energy (eV) 

R (a.u.) 
c 

Calculated 
Cl x 10 8 (cm 3 /s*c) 

Experimental 
a x 10 8 (cm 3 /sec) 

H" + 

H + 

0.7542 

24 

90 ± 34 

39 ± 2i a 

i r + 

He + 

0.7542 

26 

59 ± 20 

— 

rr + 

h 2 + 

0.7542 

25 

67 ± 23 

b 

47 ± 15 

0- + 

Na + 

1.478 

23 

20 ± 6 

21 ± 10 C 

0- + 

0 + 

1.478 

23 

22 ± 6 

27 ± 13® 

0- + 

N + 

1.478 

23 

22 ± 6 

26 ± 8 3 

0- + 

o 2 + 

1.478 

24 

19 ± 5 

10 ± 4 C 

Qr + 

n 2 + 

1.478 

23 

20 ± 6 

— 

Or + 

N0 + 

1.478 

23 

19 ± 6 

49 ± 20° 


a« 

Reference 

13, 

b. 

Reference 

20 

c. 

Reference 

14 
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For the atom-molecule systems, experimental data exist on the 
— + — *+• — 4 " 

D + H 2 , 0 + 0 2 , and 0 + NO reactions. Of these, there is satisfactory 

agreement between theory and experiment for the first two reactions with 

— + 

theory predicting larger reaction rates than observed. The O + NO 
reaction is puzzling, however, since the experimental result is more than 
a factor of 2 larger than theory. Within the framework of this theory, 
there is no satisfactory explanation of this result. 

For all the above cases, the reaction radius R is displayed in Table I. 

c 

In these reactions, the most favorable transfer of the electron occurs at 
about 25 a.u. or at AE 1 eV. The electron affinities were taken from 
the review by Berry. 21 

Molecular Negative Ion 

In systems in which the negative ion is molecular, the main difficulty 
in applying the absorbing sphere model is the determination of the correct 
electron detachment energy. In these cases, the negative ion may reside 
in an excited rotational or vibrational level or possibly in a metastable 
electronic state. Experimentally, how the negative ion is formed largely 
determines this final state. Compounding this problem is the fact that 
accepted values for electron affinities for most molecular species of 
atmospheric interest are generally not well known. Consequently, we have 
performed the calculations with the latest values for the electron affinities, 
which should be reasonably close to the real electron affinities. 

Of considerable interest to the chemistry of the upper atmosphere 
are the neutralization reactions of 0 2 . For the calculations on the 
0 2 systems, we have used a recently obtained value for the electron 
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affinity. 22 In Table II, we find both merged-beam and afterglow results 
from the 0 2 ” + 0 2 reaction. The theoretical value of 24 x 10~ 8 cm 3 /sec 
for an electron detachment energy of 0.43 eV is almost a factor of 2 below 
the merged-beam value. However, the theoretical reaction rate is in 
agreement with an afterglow experimental value of 10 x 1(T 8 cm 3 /sec. 

In this latter case, the negative ions are formed at thermal energies so 
that the electron detachment energy should be close to the electron affinity. 

In the merged-beam experiments the ions are formed in dualplasmatron discharge 
sources so we might expect excited negative ions. For example, when the 
negative ions were formed by electron bombardment in an experiment by 
Burch, Smith, and Branscomb, 23 these authors found an electron detachment 
energy of 0.15 -0.5 eV, which probably means that the ions were formed 
in an excited vibrational level. The reaction rate can then be expected 
to be larger than the calculated value, and calculations show the reaction 
rate increasing by 40$> when the electron affinity is decreased to 0.15 eV. 

It must be remembered, however, that the error in the matrix element generator, 
Eqs. (lo) and (ll), is unknown for this low value of the electron detachment 
energy. For 0 2 ~ + N 2 + , the theoretical value of 25 x 1(T 8 cm 3 /sec agrees 
well with the experimental merged-beam value. This is puzzling when 

_ 4 * 

compared with the 0 2 + 0 2 case, but it may be caused by a lack of excited 

states available for reaction in the critical R region. Several other 0,” 

c 

reactions of upper atmospheric interest are shown. 

For the N0 2 “ reactions, we have used the lower and upper limits of 
the electron affinity given by Berkowitz, Chupka, and Gutman. 24 
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ION- ION THERMAL ENERGY RECOMBINATION REACTION RATES: 


MOLECULAR NEGATIVE ION 


System 

Electron 
Detachment 
Energy (eV) 

R c (a.u.) 

Calculated 
at x 10 8 (cm 3 /sec) 

Experimental 
Ctf x 10 8 (cm 3 /sec) 

0 2 + 0 2 

0 .43 

36 

24 ± 8 

42 ± 13 a , 10 ± l b 

<V + n 2 + 

0.43 

36 

25 ± 8 

c 

16 ± 5 

0 2 ~ + 0 + 

0.43 

35 

29 ± 9 



0 2 ” + Na + 

0.43 

36 

26 ± 9 


0 2 ’ + N0 + 

0.43 

36 

24 ± 8 

— 

no 2 - + 0 2 + 

1.83 

22 

13 ± 3 

41 ± 13 S 


2. 32 

20 

12 ± 3 


N0 2 " + NO + 

1 . 83 

22 

14 ± 3 

51 ± 15®, 17.5 ± 6 d , 21 


2.32 

20 

12 ± 3 


N0 3 " + 0 2 + 

2.48 

19 

11 ± 2 

13 ± 4 f 


3.20 

18 

10 ±2 


N0 3 " + NO + 

2.48 

19 

11 ±2 

81 ± 23 f , 3.4 ± 1.2 b 


3.20 

18 

11 ±3 



a. 

b. 

c. 

d. 

e. 

f. 


Reference 13. 

Reference 16. 

Deduced from the extrapolation of data from Reference 12. 
Reference 15. 


From Reference 25 with errors estimated from the extreme spread of the data. 
Reference 20. 
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We find that the theoretical calculations here underestimate the merged- 
beam results by at least a factor of 3 for the two reactions shown. Little 
is known about the excited state populations produced in the beam experi- 
ments so we can only speculate that N0 2 ~ is possibly being formed in an 
excited electronic or vibrational state or that there is a gross deficiency 
in the theory. Another possible explanation for these large cross sections 
can be found in the work of Weiner, Peatman, and Berry 26 . These authors 
indicate that potential curve crossings at large distances, R > 100a o , 
contribute significantly and possibly dominate the cross sections. In 

our model, these curve crossings have a negligible effect on the cross 

+ 

sections. The N0 2 + NO reaction rate has also been measured at thermal 

energies by the afterglow method 15 and by plasma studies. 25 Here the 
agreement between theory and experiment is satisfactory. 

Only preliminary data are available on the "merged-beam" results for 
the N0 3 ~ reactions. For the cases available, there is agreement on the 
N0 3 - + 0 2 + reaction but an underestimation by theory in the N0 3 ~ + N0 + case. 
The afterglow results for this latter system are below the theoretical 
results and in agreement with the limitations on the absorbing-sphere 
model. In the N0 3 “ system, the electron affinity limits were again taken 
from the work of Berkowitz, Chupka, and Gutman. 24 
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CROSS SECTION ENERGY DEPENDENCE 


The energy dependence of the cross sections may be calculated 
by using Eq. (7) for the evaluation of In Fig* 1 are shown the 

predicted cross section curves for the negative ions studied here. 

The curves for each system are essentially independent of the mass of 
positive ion (within - 10 $) when the mass of the positive ion is within 
a factor of 2 to that of the negative ion. 

As expected from the previous calculations, the system with the 
lowest electron affinity, 0 2 , has the largest cross sections at a given 
energy. Likewise, the system with the largest electron affinity, N0 3 , 
has the smallest cross sections. If the negative ion lies in an excited 
level above its ground state, we would expect the observed cross sections 
to be larger than those predicted here. 

We may test the predicted cross sections with the many 0” systems 
that have been observed in the merged-beam work. This system was chosen 
because of the abundance of available data. No comparisions were made 
with the molecular negative ions since the reactant excited state populations 
are undetermined. Figure 2 draws the predicted cross section and the data 
from the five systems that have been measured. The energy dependence of 
the cross sections is very well reproduced by theory. All the cross 
sections lie well within a factor of two with theory with the exception 
of N0 + + 0* which is above the prediction throughout the energy range. 

As in the calculation of the reaction rates, the reason why the N0 + ion 
produces such large cross sections is unexplained. However, on a relative 
scale, the energy dependence of the N0 + + 0 cross sections would be well 
reproduced. 
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DISCUSSION AND CONCLUSIONS 


In the preceding sections, we have presented a model for estimating 
the upper limit cross sections for the ion-ion mutual neutralization 
reactions. The model is based on the Landau-Zener approximation, a 
semiempirical method for estimating the coupling matrix elements, and the 
assumption that there are a large number of states available for reaction. 

If there are few states available, the absorbing sphere model will definitely 

+ — 

overestimate the reaction rate, as it did in the H + H case. 

The model was set up specifically for application to molecular systems 
in which there are large numbers of neutral excited electronic and vibrational 
channels available for reaction. For systems including negative molecular 
ions, comparisons with merged-beam results were somewhat ambiguous because 
of the uncertainty in the excited state populations of the negative molecular 
ion beams. For the neutralization reaction rates that were measured at 
thermal energies by the afterglow method, the agreement with theory was . 
satisfactory. For atomic negative ion systems in which the excited state 
populations of the negative ion are not a factor, agreement with the 
merged-beam results generally was good. 

A few comments should be made concerning the significance of the 
dissociation reaction 

A + BC + -A+B+C+AE (l5) 

since this may possibly be the source of the theoretical and experimental 
differences in the molecular systems. Bates and Boyd 6 have presented the 
Landau-Zener formulas for this process, and we find that the absorbing- 
sphere model will be equally applicable to this type of reaction. The only 
difference between the theory for reaction (l) and that for reaction (15) 
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is in the energy dependence of the cross sections* At low energies, which 
is the region of interest for the reaction rates, both reactions (l) and 
(l5) will have a cross section dependence approaching E“ x . At higher 
energies, however, the cross sections for reaction (l) continue to decrease, 
but at a rate slower than E -1 , until extremely high energies of E ^ 1 keV 
are reached. The cross sections for reaction (15) at high energies, however, 
approach a constant value- The cross section energy dependence at high 
energies will be determined primarily by the neutral molecular states of 
the positive ion* We may predict that reaction (15) will be important if 
the neutral molecular states of the positive ion in the potential energy 
region of the ground state positive ion minus the electron detachment 
energy of the negative ion minus approximately 1 eV lead to dissociation* 

So far, the energy dependence of the experimental data on the molecular 
ions tends to indicate that the dissociation reaction is not a dominate 
process for those systems measured. 

With regard to the vibrational change of state on reaction, the 
Franck-Condon factors, Eq. (14) will dictate that at large distances the 
transition probability for the transfer of the electron will be greatest 
for states with the greatest amount of vibrational overlap. Generally this 
means that for reactions at large distances, which determine the magnitude 
of the cross sections the most, the conditions are favorable for little 
change in the vibrational state. At smaller intermolecular separations, 
where the electronic matrix coupling element alone is so large that it 
precludes reactions, transitions to states where there is little vibrational 
overlap will allow a significant amount of reaction. Generally speaking, 
we would expect a greater change in vibrational energy to occur for 
reactions at small impact parameters. 

Probably the most significant factor to note in this work is the 
large reaction rate dependence on the electron detachment energy. If the 
energy for the electron detachment is decreased because of populated 
excited states of the negative ion, the cross sections increase. If the 
excitation is to the upper levels, this increase becomes significant 
For upper atmosphere chemistry, this fact necessitates knowing the excited 
state populations of the reactants for the ion-ion recombination reaction. 
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FIGURE CAPTIONS 


1* Cross section energy dependence predicted by the "absorbing-sphere" 
model. For each system, the cross section curve is accurate to 
± 10$ for mass ratio variations of less than 2. The electron detachment 
energies used were: 0 2 - 0.43 eV, H - 0.754 eV, 0-1.48 eV, N0 2 - 2.1 eV, 
and N0 3 - 3.0 eV. 

2. Solid line: cross section curve predicted by the "absorbing-sphere" 
model for 0 systems. The electron detachment energy used for 0 was 
1.48 eV. The data are from "merged-beam" results; references 9, 12, and 
14. 
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Abstract 

The ionization (Penning and associative) total cross sections for 
Ne*(2p 5 3s 3 P) + Ar are calculated semiclassically in the energy range of 
0.001 to 1000 eV and compared with recent experimental results of Micha, 
Tang, and Muschlitz. 
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The recent experimental data of Micha, Tang, and Muschlitz 1 on the 
low energy (0.01 to 0,1 eV) dependence of the ionization total cross sections 
for the sum of the Penning ionization process 

Ne*( 3 P 2 , 3 P 0 ) + Ar - Ne + Ar + + e~ (l) 

and of the associative ionization process 

Ne*( 3 P 2 , 3 P 0 ) + Ar NeAr* + e (2) 

display an unexpected energy dependence that has not been explained by 
theoretical calculations. These authors find cross sections that first 
decrease as the energy increases, then pass through a minimum, and later 
increase at the highest energies measured. Calculations by these authors 
using a Lennard-Jones, L-J, (l2,6) interaction potential and a short range 
coupling width T(r) were able to reproduce the low energy, monotonically 
decreasing behavior of the cross sections but did not follow the change in 
slope. Further experimental work by Tang, Marcus, and Muschlitz 2 on 
Ne* + Kr and Ne* + Xe show the same type of behavior. On another noble gas 
system — He* (2 3 S and 2*s) + Ar — the thermal energy cross section 3 * 4 is 
comparable or smaller than the high energy results at 100 to 1000 eV. 5 It 
therefore appears that the shape of the energy dependence of the Ne* + Ar 
cross sections could be common to many of the noble gas ionization reactions. 

Calculations were performed using the semiclassical formulas found in 
a paper by Miller 6 . Here, the ionization cross section is 

05 

Q = 2tt<*J db'b'P(b) (3) 

o 

where the ionization probability is given by 
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The radial velocity is calculated from 

v b (H' = V C [ i - v 0 ;rV'e - b a /n a y (s', 

"where v c is the incident velocity and b is the impact parameter. In 
Hq. '3 \ ;v is the statistical weight factor that is unity for this system. 

"-"he coupling to the continuum is provided by the width r(R), which in this 
case ve have assumed to be of exponential form. This form is similar to the 
shapes of calculated bv Hiller, Slocomb, and Schaefer'; Bell 8 ; and 

Fuji! et al? for the He*' + II system and is similar to that deduced from 
He*’ + l T a experimental data by Hotop and Niehaus 10 . The interaction potential 
V c ( R ) is that of the rnetastable atom in the field of the neutral atom. Severn 
different forms for the interaction potential were used while we restricted 
the coupling width to the exponential form 


r(R) = A exp ( -BR) . (6) 

V/c first hoped that the data could be reproduced by a potential with 
two minima, arising from, a long range van dor Waals well, a barrier at 
intermediate distances with a magnitude equal to the energy at which the 
experimental cross sections change slope, and a strong chemical well at small 
separations. In Figure 1, the short dashed lines show such a potential that 
v/as formed to reproduce approximately the L-J potential of Micha et al. 1 
**rom R ^7.5 to 10 a D ^ Its functional form is given by 


V 


o v 


(R) 




+ a 



(7) 
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and the parameters used were = 0.15 a*u. , c = 3.5 a.u. , a = 2.5 x 1CT 3 a.u., 

b - 1.0 ii~ l , and c = 6.75 a 0 . A second calculation was made using the L-J 
potential parameters of Ref. 1 that were obtained from the Slater-Kirkwood 
value l'or the long range dipole-induced dipole term and the composition rules 
for <;. T . The parameters are e = 1.76 x 10“ 4 a.u., o = 8.16 a D , and a = b = c - 0. 
This L-J potential is shown by the solid line. 

In Figure 2 are displayed the experimental data and the results for the 
L-J potential, solid line, and for the potential with two minima, dashed 
line. The coupling width, Eq. (6), had the parameters A =0.01 a.u. and 
B = 1.0 a” 1 . For collision energies greater than the barrier height, 0.04 eV, 
we see the cross sections rising dramatically when compared with the L-J 
results. Since the experiments did not show such a rapid rise, the width 
was varied to try to obtain better agreement. We reduced the width for 
separations less than the potential barrier and set it equal to its value at 
R = 6.9 a 0 for R < 6.9 a 0 . The results are given by the short-long dashed 
line, but again there is a sharp rise at E ~ 0. 04 eV that is not borne out 
by the data. 

Another possibility to explain the data is that t(r) should be increasing 
faster at small R so that as the particles penetrate to smaller separations 
at the higher energies, the transition probability will increase. Rather 
titan varying F(r) randomly, we set A = 1.0 and chose B such that it best 
reproduced the low energy data around v = 500 m/sec. When the L-J potential 
oi Ref. 1 was used, we obtained the results shown in Figure 3 by the solid 
line. The physical characteristics of the cross sections are well reproduced, 
and the value of B obtained was 1.53 a 0 “ 1 , close to the exponential character 
of various calculated widths 7- " 10 that vary between 1 and 2. 

Since the L-J potential possesses a much too repulsive core for the 
higher energy work, the calculations were also performed using a Morse potential 
with the same minimum position and magnitude as the L-J potential. It is given 

e -2-(l-R/Re) _ 2e - 3(l-R/Re) ") (8) 
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where ~ = 1.76 x 10“ 1 a.u. , R =9.16 a n , and - =6.6. The constant - 

e u 

was chosen to reproduce approximately the repulsive potential of Abrahamson 1 1 
for Xa Ar in the region of 1 to 3.5 a 0 where V q (r) decreased from approxi- 
mately 1000 to 10 eV. The L-J potential in this region is 2 to 3 orders 
of magnitude too large. The Na 4- Ar potential was used since it would be 
closest to reproducing the electron distribution of Ne* + Ar. 

V.’ith the Horse potential, the low energy behavior is similar to that 
of the L~J potential but now at higher energies the cross sections are 

c n 

much larger, ~ 10 A , and similar in magnitude to the He* -f Ar measurements. 5 
Since the Horse potential has a more realistic repulsive core, we would 
expect that calculations made using it would be more accurate at high 
energies than the L-J results. Further experimental measurements will be 
needed to test the energy dependence and the magnitude of the high energy 
cross sections predicted by the Morse potential calculations. Also, since 
the width has been normalized to the Ne* + Ar low energy cross sections 
that arc uncertain by 50^, it would be unrealistic to expect that the 
magnitude of the theoretical calculations at higher energies will be more 
accurate than a factor of 2. 

Most important, however, is the fact that only one parameter in the 
coupling width was varied to obtain the correct normalization to the 
experimental data at thermal energies. Since the interaction potential 
between the metastable and neutral atom can be derived from other sources, 
estimates of the energy dependence of the ionization cross sections can be 
made by using the available thermal energy results. 
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FIGURE CAPTIONS 


Potentials used in the calculations: L-J (12,6) - solid line, 

Morse - long dashed line, and the double minima potential - short 
dashed line. 

Ionization total cross sections calculated using the L-J potential 
and ~(R; = 0. 01 * exp (-R) is given by the solid line. For the double 
minima V(r) and the same j~(R), the results are shown by the dashed line, 
and for the double minima V(r) and the truncated F(r) the results are 
displayed by the short-long dash line. The experimental data of Ref. 2 
are shown by the solid circles. 

Ionization total cross sections calculated with F(r) = exp (-1.53R) 
and the L-J potential, solid line, and with the Morse potnetial, 
dashed line. The experimental data of Ref. 1 are shown by the solid 
circles. 
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